COUNTING MONOMIALS 



MORDECHAI KATZMAN 

Abstract. This paper presents two enumeration techniques based on Hilbert functions. The 
paper illustrates these techniques by solving two chessboard problems. 



1. Introduction and preliminaries. 

The purpose of this note is to illustrate two powerful enumeration techniques based on computa- 
tional Commutative Algebra methods. 

By way of illustration I chose to apply these methods to the following two elementary problems: 

(1) Consider a n X n chessboard. What is the maximal number of unattacked squares in the 
board after placing on it k queens? More generally, in how many ways can we place k queens 
on a chess board to obtain exactly u unattacked squares? 

(2) Consider an infinite chessboard. How many squares can a knight reach in d moves? How 
many squares can be reached in d moves and no less? 

Although these problems are phrased in the language of chess, they are specific instances of more 
general graph-theoretical problems. The enumeration techniques presented here answer these more 
general problems. 

At the heart of the methods presented in this paper are the notions of graded modules and their 
Hilbert functions. In essence, we will reduce each of the problems above to a problem about the 
enumeration of sets of monomials, and this enumeration will be achieved using Hilbert functions. 

While the application of Hilbert functions to the problems presented in this paper is new, the 
use of Hilbert functions in combinatorics is not. The solution of some simple enumeration problems 
using Hilbert functions, such as finding the independence number of a graph, has long been part of 
the folklore of computational commutative algebra experts. An early and striking example of the 
use Hilbert functions in combinatorics is Richard P. Stanley's work on magic squares (I refer the 
reader to [Hj for an accessible and thoroughly enjoyable account of this work.) 

We now review graded modules and Hilbert functions. Throughout this paper, all rings are 
commutative and with 1; K will always denote a field. 
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A if-algebra R is N^-graded if we can write 

R = ^a, 

a direct sum of abelian groups, and the direct summands satisfy 

RaRb C i? a +b 

for all a, b <E N N . Henceforth we shall also impose the condition R = if, which implies that each R a 
is a if- vector space and that, if R is a finitely generated if-algebra, each R a is a finite dimensional 
if -vector space. For each a G N w we shall refer to the elements of R a as being homogeneous of 
degree a. 

A fundamental example of such a graded if-algebra is the ring of polynomials R = if [xi, . . . , x n ). 
We can endow R with different graded structures. We are all familiar with the N-grading 

i? = 0i? a 

in which each R a consists of the homogeneous polynomials of degree a. We can define another 
grading as follows: let di,. ..,d„ G N w and define the degree of a monomial x" 1 . . . x" n to be 
aidi + . . . a n d n . We can now write 

R = Ra, 

aGN™ 

where each R a is the if- vector space spanned by all monomials of degree a G N^. 

Let R be a N^-graded if -algebra. An i?-module M is graded if it has a N w -grading compatible 
with that of R, i.e., if we can write 

M = M B) 

a direct sum of abelian groups, and the direct summands satisfy 

R a M b C M a+h 

for all a,bGN w . 

If R is a polynomial ring as in the examples above and I C R is a homogeneous ideal, i.e., an 
ideal generated by homogeneous elements, then R/I has a natural structure of a graded i?-modulc. 

Let R be a N^-graded if -algebra and let M be a graded i?-module. We define the Hilbert function 
HF A f of M to be the function HF A f : N w -> N defined by HF M (a) = dim K M a . The Hilbert series 
HSm(£ij • ■ • , ijv) of M is the generating function of the Hilbert function, i.e., 

HS M (h,...,t N )= HF m (a)t? 1 • • • *7 • 
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If R is a polynomial ring as in the examples above with its familiar N-grading, and if we view 
R as a graded i?-module, then HF r(o) is just the number of monomials of degree a in n variables, 
i.e., HF^(a) = and HSjj(i) = 1/(1 - t) n . If we were to assign degrees di, . . . , d n e to 

x\, . . . , x n we would obtain 

HS_r(£i, . . . , ijv) 



n n 1 _ + dil +d iJV ' 
»=1 1 l l • ■ ■ l Ar 

Take i? to be a polynomial ring with its familiar N-grading, let J C -R be a homogeneous ideal and 
write S — R/I. One can show that HFs(a) is of polynomial type, i.e., it agrees with a polynomial, 
the Hilbcrt polynomial HPs(a) of 5, for all a>0. The degree of HP5 is one less than the Krull 
dimension of S. Also, one can write 

TTO M 

where P(t) is a polynomial which does not vanish at t = 1 and d is the Krull dimension of S. 

2. Unattacked squares 

We now consider the first question mentioned in the introduction. We naturally identify the 
squares of the n x n chessboard with pairs (i, j) where 1 < i,j < n. 

We fix n, the size of the board. Let K be any field and define R to be the polynomial ring in 2n 2 
variables 

R = K[xu, . . . ,x nn ,yu, ■ ■ ■ ,y nn ]- 

We assign degree (1, 0) to all the x variables and degree (0, 1) to all the y variables. 

Roughly, the x variables will correspond to squares in our n x n chessboard which are occupied 
by queens while the y variables will correspond to unattacked squares on the board. 

We define / to be the ideal of R generated by the squares of all variables together with 

{xijyim I a queen can move from square (i, j) to square (l,m)} . 

Notice that /, as any other ideal generated by monomials, is homogeneous with respect to the 
N 2 -grading of R. 

For any k > define 

H{k) = max{/i G N | dim K (R/I) {k ^ > 0}. 

Proposition 2.1. fi(k) is the maximal number of squares on the nxn chessboard which can remain 
unattacked after placing on it k queens. 
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Proof. Consider any monomial M = x^y' 3 in R whose image in R/I is not zero. Since / contains 
the squares of all the variables, M must be square-free and we may write 

where all the variables in this expression are distinct. We next observe that for any 1 < £ < A and 
1 < C < v, a queen cannot move from square to square (l^,m^), otherwise, j £ y; c ,m c would 

be one of the generators of / and M would be zero modulo /. We showed that every monomial of 
degree (A, fi) whose image in R/I is not zero corresponds to a configuration on the chessboard where 
the squares ■ ■ ■ , are occupied by queens and the squares {l\,m\), . . . , (l v ,m v ) arc not 

attacked by any of these queens. 

It is easy to see that the converse is also true and so we have established a bijection between the 
configurations of A queens and v unattacked squares and the set of monomials of degree (A, v) which 
are not zero modulo /. 

Notice that all the graded components (R/I)^ x ^ are spanned as if- vector spaces by monomials 
of degree (A, v), and that a basis for (R/I)^ x ^ is given by the set of all such monomials whose 
images in R/I are not zero. So now we can see that the condition 

dim K (R/I) M > 0, dim K (iJ//) (fci „ +1) = 

can be translated using the bijection established above to the statement that it is possible to place 
k queens on the chessboard so that one can find [i unattacked squares but not /i + 1 unattacked 
squares. □ 

We now address the more general question: in how many ways can we place k queens on 

a chessboard to obtain exactly u unattacked squares? 

Proposition 2.2. For any < u < fi(k) 

n(k) , . 

$(M) = HF fl/J (M)- E LP^' 

v=u+l ^ ' 

Proof. We proceed to prove this by reverse induction of u. When u = fi(k) the equality $(fc, /J,(u)) = 
HF R /[(k, n(u)) follows easily from the discussion in the proof of the previous proposition. 

Pick now any < u < n(k). HF ~ R /i{k,u) is the number of ways one can choose the position of 
k queens and u squares unattacked by these queens. For each such choice, one can extend the set 
of u unattacked squares to a maximal set of v unattacked squares by the same k queens. To obtain 
&(k, u) we need to count only those choices for which u = v or, equivalently, we need to subtract 
from HF R /j(k,u) the number of configurations which which extend to a maximal one with v > u 
unattacked squares. The induction hypothesis implies that there are exactly configurations 
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with k queens and a maximal set of v unattacked squares, and each one of these produces 
configurations with k queens and u unattacked squares which can be extended to a maximal set of 
v unattacked squares. Subtracting all these, we get the desired result. □ 

Table 1 lists the values of <E>(fc,w) when n = 8 for 3 < k < 43 and 1 < u < 25 (blank entries 
are zero.) For example, the table shows that /i(8) = 11 and that $(8, /i(8)) ~ 48, which means 
that the largest number of unattacked squares one can have when 8 queens are placed on a regular 
chessboard is 11, and that there are 48 such configurations. This is the answer to a question 
originally published by W. W. Rouse Ball in 1896 |2| (see also chapter 34 in 3 .) This calculation 
was produced by FreeSquares, a CH — h program which can be found in [5]. (There are several widely 
used computer packages which can compute multi-graded Hilbcrt series, but unfortunately they are 
not very efficient.) 

The method introduced in this section generalizes naturally to deal with graph-theoretical prob- 
lems which we now describe. Let G be a finite graph. If U and W are disjoint sets of vertices of G 
we say that U and W are independent if there is no edge connecting a vertex in W with a vertex 
in U. For a given k what is the maximal size of a set of vertices which is independent of a set of k 
vertices? In how many ways can one choose independent U and W with given size? 

Let {v i, . . . , vn} be the vertices of G. One obtains the solution to this more general problem by 
replacing the ring R with K[x±, . . . , xn, y%, ■ ■ ■ , Un] and the ideal / above with the ideal generated 
by the squares of all the variables and 

{xiy 3 | (vi, Vj) is an edge in G} . 

3. Knight moves in an infinite chessboard. 

We now consider the second set of questions mentioned in the introduction: How many squares 
can a knight in an infinite chessboard reach in d moves? How many squares can be reached in d 
moves and no less moves? We will denote the first number with f(d) and the second with g(d). 

The implementation of the results in this section relics on Grobner bases techniques- the reader 
may want to consult for an introduction to Grobner bases. However, to appreciate the general 
ideas behind the approach of this section no knowledge of Grobner bases is needed. 

We again let K be any field and let R be the if-subalgebra of K [x\, X2, , x^ 1 } generated by 

M_ f 2 2 -12-2 -2 2-1-1-2 -2 -1\ 

— I 1 3> 2 5 *^ \^- ,! 2, 7 *^ 2. *^ 2 ' ■2'2;*^1"^'2 ) *^ 1 5*^1 2 ) 1 *^2 f* 
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Table 1: Values of <&(&,«): rows correspond to values of k 
while columns range over values of u. 
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The first step towards the solution of this problem is to realize that f(d) is the cardinality of 
M d :— {ai . . . ad\ai, . . . , ad € M} while g(d) is the number of elements in M d but not in any M l for 
i < d. 

We can produce a presentation for R by mapping a polynomial ring S = K[y\, . .. to i? by 
Hi — > rrii where rrii is the ith element of M. We denote this mapping with if>. Notice that the 
restriction of iff to the set of degree-d monomials in S gives a surjection onto the elements of M d . 

Let k be the kernel of the map above. This kernel can be computed effectively using Grobner 
bases techniques as follows: let I be the ideal of k[u, xi,x 2 ,yi, ■ ■ ■ , ys] generated by 

{uxix 2 - l,yi - Xj_xl, y 2 - x\xi,y z XY - x\,y±x\ - x 2 , 
y 5 x 2 - xi,y 6 x 2 - x\, yix x x\ - 1, y%x\x 2 - 1} 

and fix an elimination order where u, x±,X2 are the largest variables. Then n is generated by the 
elements of a Grobner basis for / which do not contain the variables u, x\, x 2 (cf. chapter 1 of 7 .) 
Recall also that k is a binomial ideal. 

Notice that the ring R is not very interesting: it is in fact identical to K \x\, x^ , x 2 , xl^ ] (here 
is a chess proof: x\ G R because a knight can move one square to the right in three moves. By 
symmetry also x± , x 2 , x^ 1 £ R-) However, S/k is far more interesting for reasons explained below. 

Since the restriction of iff to the set of degree-d monomials in S is a surjection onto M d , to find 
f(d) we need to find the size of a maximal set of degree-d monomials in S which are distinct modulo 
k. Two such monomials y a and y@ are distinct modulo n if and only if y Q — y^ 3 is not in the largest 
homogeneous sub-ideal H of n. It is easy to compute H: the elements of H are the elements of 
the homogenization of k with respect to a new variable, say t, which do not involve t, thus we can 
compute H by homogenizing a Grobner basis for K using a graded lexicographic order (cf. exercise 
1.6.19 in pQ) and eliminating the variable t. We notice that this Grobner basis can be chosen to 
consist of binomials, and so H is also a binomial ideal. 

So we have reduced the problem of computing f(d) to the problem of finding the size of a maximal 
set of degree-d monomials in S which are distinct modulo H. Fix any term ordering in S and let TL 
be a Grobner basis for H consisting of binomials. Now for any two monomials y Q > y^ of the same 
degree, y a = y 13 modulo H if and only if y a reduces to y^ with respect to TL. Since each reduction 
of a monomial with respect to TL produces a new monomial (of same degree) , to produce a maximal 
set of degree-<i monomials in S which are distinct modulo H we may pick all monomials of degree 
d which are non-zero modulo m(H), i.e., 



f(d) = dim K (S/m(H)) d = dim K (S/H) d = UF S/H (d) 
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where the second equality is a celebrated theorem proved by F. S. Macaulay in [Hj. 
An easy computation with Macaulay2 (|3j) shows that 

uo u . 1 + 5* + 12i 2 - 8t 4 + At 5 
tib S /H{t) (1 -t) 3 

and that the Hilbert polynomial of S/H is 1 + Ad + 7d 2 . Since 



HS s /jr(t) - ^(1 + Ad + 7d 2 )t d = -At 2 - At 



we obtain 



f(d) = { 



1 d = 

8 d = 1 

33 d = 2 

1 + Ad + 7d 2 d > 3 



We now proceed to compute g(d). We again fix a monomial ordering in S which refines the total 
degree ordering. List all the monomials in S in ascending order, and let B be the set of all degree-d 
monomials in S which are not congruent modulo k to a monomial appearing earlier in the list. We 
now show that g(d) = #£>. 

If for two distinct degree-d monomials y a > y 13 we have f (y a ) = ^(y' 3 ) then y Q — y 13 £ k 
contradicting the choice of B. Hence the restriction of $ to B is injective. Similarly, if for some 
degree-d monomial y a there exist a monomial y 13 of degree i < d so that \&(y a ) = ^(y ) then 
y<* _ yP £ K anc i since y a > y 13 we get a contradiction to the choice of B. Hence the restriction of 
^ to B is a surjection onto M d \ Ui^M 1 . 

Using the fact that K has a Grobner basis generated by binomials we may deduce that B is the 
set of all monomials which are not in in k and so 

g(d) = dim K (S/m(K)) d = HF (s/in(K)) (d). 

Another straightforward computation with Macaulay2 shows that 



1 + 6t + 17t 2 + 12t 3 - 8t 4 - At 5 + At 6 



HS (s/in(K)) (t) _ 



and that the Hilbert polynomial of S/ in(«;) is 28d — 20. Since 

oo 

HS( S/ln(K) )(t) - J^(28d - 2Q)t d = At* + At 3 - At 2 + 21 



i=0 
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we obtain 



9(d) 



1 


d = 


8 


d= 1 


32 


rf = 2 


68 


rf = 3 


96 


rf = 4 


28d - 20 


d > 5 



The methods of this section also generalize in a natural way. Let 



W = 



( W\\ \ 



\ Wlm J 



( WN1 \ 



\ W Nm J ^ 



c 



be a finite set and consider an infinite directed graph G whose vertex set is Z m and for any u,v € Z m , 



(u, v) is a directed edge if and only if v — u S W. 

By replacing R above and its presentation S — * R with the presentation 



K[ yi ,...,y N }^K[x^ C-V..,^ 1 



which maps yi to x™" m for all 1 < i < N, we can, by following exactly the same procedures 

as before, produce closed formulas for the functions f(d) which count how many endpoints all length 
d paths starting at a fix vertex have, and closed formulas for the functions g(d) which count how 
many vertices are at a distance of d from a fixed vertex. 

Theorem 3.1. For any directed graph G as above, there exist polynomials P(d) and Q(d) so that 
f(d) = P{d) and g{d) = Q{d) for all d > 0. 

Proof. This is an immediate consequence of the fact that Hilbert functions are of polynomial type. 

□ 



Appendix: A Macaulay2 implementation. 

All the methods in this paper are easy to implement with existing computer systems. As an 
example aimed to tempt the reader to experiment with these systems we present a Macaulay2 
program for the solution of the enumeration problem in the previous section: 

R=ZZ/101 [u,a,b,y_{l>. .y_{8},Monomial0rder=>Lex] ; 
I={u*a*b-l_R,y_{l}-a*b~2,y_{2}-a~2*b,y_{3}*a-b~2,y_{4}*a~2-b, 
y_{5}*b"2-a,y_{6}*b-a~2,y_{7}*a*b"2-l_R,y_{8}*a~2*b-l_R}; 
G=gens gb ideal I; 
J=selectInSubring(3,G) ; 

Sl=ZZ/101[y_{l>. .y_{8},t] ; 
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J=substitute ( J , SI) ; 
H0=homogenize (gens gb J,t); 

S2=ZZ/101 [t,y_{l>. .y_{8},Monomial0rder=>Lex] ; 
H0=substitute(H0,S2) ; 
G=gens gb ideal HO; 
H=selectInSubring ( 1 , G) ; 

S=ZZ/101[y_{l}. .y_{8>] ; 
J=substitute ( J , S) ; 
H=substitute(H,S) ; 
print (hilbertSeries coker J); 

print (hilbertPolynomial (coker J, Projective=>f alse) ) ; 
print (hilbertSeries coker H) ; 

print (hilbertPolynomial (coker H, Projective=>f alse) ) ; 

This produces the following output: 
6 5 4 3 2 
4$T -4$T -8$T +12$T +17$T +6$T+1 



2 

(-$T+1) 

28$i-20 

5 4 2 
4$T -8$T +12$T +5$T+1 



3 

(-$T+1) 

2 

7$i +4$i+l 
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